Weak localization and conductance fluctuations in a quantum dot with parallel 
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In the presence of both spin-orbit scattering and a magnetic field the conductance of a chaotic 
GaAs quantum dot displays quite a rich behavior. Using a Hamiltonian derived by Aleiner and 
Farko [Phys. Rev. Lett. 87, 256801 (2001)] we calculate the weak localization correction and 
the covariance of the conductance, as a function of parallel and perpendicular magnetic field and 
spin-orbit coupling strength. We also show how the combination of an in-plane magnetic field and 
spin-orbit scattering gives rise to a component to the magnetoconductance that is anti-symmetric 
with respect to reversal of the perpendicular component of the magnetic field and how spin-orbit 
scattering leads to a "magnetic-field echo" in the conductance autocorrelation function. Our results 
can be used for a measurement of the Dresselhaus and Bychkov-Rashba spin-orbit scattering lengths 
in a GaAs/GaAlAs heterostructure. 
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I. INTRODUCTION 

A two-dimensional (2D) electron gas offers ample op- 
portunity for manipulation of the electron's spin and or- 
bital state. In the absence of spin-orbit coupling, such 
manipulation can be performed separately on the elec- 
tron spin and orbital degree of freedom, by means of mag- 
netic fields of various orientations and electrostatic gates. 
On one hand, a magnetic field parallel to the plane of the 
2D gas lifts the spin degeneracy and allows for a measure- 
ment that distinguishes between the transport properties 
of "up" and "down" spins: As long as the confining po- 
tential that determines the 2D electron gas is sharp, it 
does not significantly affect the electron's orbital degrees 
of freedom. On the other hand, gate voltages and a weak 
magnetic field perpendicular to the quantum well act on 
the orbital degrees of freedom and do not couple signifi- 
cantly to the electron's spin.^ 

Spin-orbit scattering couples the spin and orbital de- 
grees of freedom. Conventionally, the main signature of 
spin-orbit coupling in the 2D electron gas in GaAs het- 
erostructures is the observation of weak anti-localization, 
a small positive quantum correction to the conductiv- 
ity that is suppressed by a magnetic field. ^ Weak anti- 
localization is the counterpart of weak localization, a neg- 
ative correction to the conductivity at zero magnetic field 
caused by the constructive interference of back-scattered 
electron waves in a phase-coherent disordered conductor 
in the absence of spin-orbit scattering.'^ Only recently, 
spin-dependent phase coherent transport was considered 
for a quantum dot, a small island of the electron gas con- 
fined by gates'*. Motivated by an experiment by Folk et 
al.,^ theoretical works by Halperin et al.^ and by two 



of the authors,^ have shown that the combination of 
spin-orbit scattering and a spin degeneracy lifting par- 
allel magnetic field leads to a remarkably rich structure 
of quantum interference phenomena in III-V semiconduc- 
tor quantum dots, that far surpasses a simple interpola- 
tion between weak localization and weak antilocalization 
physics. 

Not only does the crossover regime between weak lo- 
calization and anti-localization represent an interesting 
issue of its own, understanding features imposed upon 
quantum transport by spin-orbit coupling is needed for 
a correct extraction of the decoherence time from exper- 
imental data. For electrons in a small dot with weak 
but finite spin-orbit coupling, the the crossover of weak 
localization to anti-localization manifests itself as a sup- 
pression of the zero-temperature value of localization cor- 
rection to the dot conductance. This suppression of the 
weak localization correction may be misinterpreted as a 
saturation of the decoherence time at low tempera- 
tures. Without a priori knowledge of the spin-orbit cou- 
pling constants for a given semiconductor dot, one must 
use the features of spin-orbit coupling induced interfer- 
ence effects in order to identify the origin of what can be 
mistaken for a suppression of the interference part of con- 
ductance, when the observed weak localization correction 
is less than the prediction of the quantum transport the- 
ory for to the crossover between orthogonal and unitary 
symmetry classes. 

In this publication, we present a detailed quantitative 
analysis of the influence of the coupling between elec- 
tron spin and orbital motion in a two-dimensional semi- 
conductor on the interference corrections to the conduc- 
tance of a quantum dot. We consider the effects of en- 
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hanced/supprcsscd back-scattering and the variance and 
correlation properties of universal conductance fluctu- 
ations. The calculations, which are performed in the 
framework of the random scattering matrix approach of 
random matrix theory, are described in Sections III and 
IV. In Section III, we present the analysis of the average 
conductance of a chaotic dot. Section IV is devoted to 
conductance fluctuations. For technical reasons, we limit 
our attention to the case of a large number of channels, 
A^i and N2 in the point contact connecting the dot to 
the bulk, though all qualitative features we discover for 
iV ^ 1 would persist in the case of TV ^ 1. These two 
technical sections are preceded by a qualitative discus- 
sion (Sec. II) and followed by an analysis of the effect 
of a spatially non-uniform spin-orbit coupling strength 
(Sec. V). 



II. INTERPLAY BETWEEN SO COUPLING 
AND ZEEMAN SPLITTING IN QUANTUM 
DOTS 

In GaAs/AlGaAs hcterostructures, spin-orbit coupling 
owes its existence to the asymmetry of the potential con- 
fining the 2DEG and the lack of inversion symmetry in 
the crystal structure, leading to the Bychkov-Rashba and 
Dresselhaus terms in the Hamiltonian. In the theory we 
present below, we consider electrons in a hcterostructure 
or quantum well lying in the (001) crystallographic plane 
of a zinc-blend type semiconductor and choose coordi- 
nates Xi,X2 along crystallographic directions ei — [110] 
and §2 = [110] (coordinate is perpendicular to the 
plane of the two-dimensional electron gas) . The effective 
two-dimensional Hamiltonian of electrons in a quantum 
dot takes the form 
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+ V{r) + -fXBgB ■ a. 
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where B is the magnetic field, A = 53(63 x r)/2c is 

the vector potential corresponding to the component of 
the magnetic field perpendicular to the plane of the two- 
dimensional electron gas, and a = (<ti, (T2, cts) the vec- 
tor of Pauli matrices. The potential V both confines 
the electrons to the quantum dot and describes clastic 
scattering from non-magnetic impurities in the dot. The 
two length scales, Ai and A2 are associated with spin- 
orbit coupling for an electron moving along the principal 
crystallographic directions §1 and §2,^ and characterize 
the length at which spin of an initially polarized electron 
would precess with an angle 27r. 

For a dot with homogeneous electron density and, 
therefore, parameters of confining potential, spin-orbit 
coupling parameters Ai and A2 are independent of the 
position. In that case, the spin-orbit coupling takes the 
form of a spin-dependent "vector potential", which is 



non-abelian since the Pauli matrices do not commute 
with each other. If Ai and A2 arc large compared 
to the dot size L, the non-commutativity of this spin- 
dependent vector potential involves higher powers of 
i/Ai,2, so that the spin-orbit couphng can be "gauged 
out" to leading order in L/Ai,2. This may be achieved 
through the unitary transformation ^'(r) = U{r)il){r), 
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(2) 
(3) 



which performs a position-dependent rotation of the spin 
out of the plane of the two-dimensional electron gas to the 
locally adjusted frame.'' Below, we assume that i? <C 1, 
and use i? as a small parameter to derive the form of 
the Hamiltonian H = WHU in the locally rotated spin- 
frame. 
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Here we introduced the spin-dependent vector potential 
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a_L = 



4A1A2 



[63 X rj 



that has the same form as a magnetic field ±6063/2 A1A2 

with opposite directions for electrons with spins "up" and 
"down" in the new spin-frame, and we also abbreviated 

1 f Xiai X2a2\ , 
+ [63 xr] 



6A1A2 \ A 

/ BiXi B2X2\ 
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The Hamiltonian (4) describes electrons in a rotated 
spin frame. It contains all relevant terms that lift the high 
degree of symmetry of a system with uncoupled orbital 
and spin degrees of freedom, to leading order in the small 
parameter L/Ai^2- In these equations, we omitted sub- 
leading terms of higher order in L/Ai,2 that do not affect 
the symmetry of the Hamiltonian. 

As long as the rate at which electrons escape from the 
quantum dot into the leads is much smaller than the 
Thouless energy i?Th, transport properties of the quan- 
tum dot can be calculated using random matrix theory. -'^^ 
The use of random matrix theory requires an analysis of 
the symmetries of the scattering matrix, which are set by 
the relative magnitudes of characteristic energy scales for 
each of the terms in the Hamiltonian (4) and the escape 
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rate. These energy scales are 
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for the orbital contribution of the magnetic field, the 
spin-orbit terms a|| and a±, and the Zeeman coupling 
terms h'^ and h"^, respectively. For weak uniform spin- 
orbit coupling in a small dot such that -L/Ai,2 1, one 
has the strong inequalities 



sl < £^. 



(7) 



Here i?Th is the Thouless energy, which is the largest 
energy scale in the problem, = Ai A2, and k, k' , and k" 
are coefficients of order unity. The coefficients k' and k" 
depend on the sample geometry and on the ratio A1/A2 
of the spin-orbit lengths. The coefficient k" also depends 
on the direction of the magnetic field. The coefficient k 
depends on the sample geometry only; it appears both for 
the orbital contribution of the magnetic field and for the 
spin-orbit term a± because both terms have the same 
spatial dependence. For a circular quantum dot with 
radius L, mean free path / ^ L, Fermi velocity vp, and 
Ai = A2 = A, one has ETh = hvpl/^L"^, n = 2, k' — 1/3 
and k" « 0.292, see App. A. 

If the spin-orbit coupling is non-uniform, i.e., when 
the spin-orbit lengths Ai and A2 depend on position, 
the "spin-dependent vector potential" in Eq. (1) can no 
longer be gauged away to leading order in L/Xi^2- A 
spatially non-uniform spin-orbit coupling can be created 
intentionally, with the help of a metal gate parallel to the 
two-dimensional electron gas that changes the asymme- 
try of the quantum well.^^'^^'^^, or arise as a by-product 
of the confining gates. Generically, the non-uniformity 
gives rise to a term in the Hamiltonian of the same sym- 
metry as a|| in Eq. (4).^^ Hence, the main effect of non- 
uniformities is to increase the corresponding energy scale 
£|l° by an amount of order i?Th(i/Afl)^, where 1/Af is 
a measure of the fluctuations of the spin-orbit coupling 
Xi2-^^ For very small quantum dots, this increase of £||° 
may eventually reverse the first inequality in Eq. (7). 

It is the existence of two energy scales each to describe 
the strength of the spin-orbit and Zeeman terms in the 
Hamiltonian that leads to the rich parameter dependence 
of the conductance distribution of a chaotic quantum dot. 
When all relevant energy scales are either much larger or 
much smaller than the escape rate, the scattering ma- 
trix has a well defined symmetry, and the conductance 
distribution can be found using symmetry considerations 
alone. ^ It is the goal of this paper to address the gen- 
eral parameter regime, where an interpolation between 
the various symmetry classes is called for. Our results 



are important for a quantitative analysis of recent exper- 
iments by Zumbiihl et al.,^^ and for the development of 
techniques for experimental determination of spin-orbit 
coupling parameters. Looking into various regimes, we 
show how to extract such information from (a) weak lo- 
calization measurements in the presence of an in-plane 
magnetic field; (b) an analysis of how the in-plane mag- 
netic field gives rise to a component of the magneto- 
conductance that is anti-symmetric in the perpendicular 
magnetic field; and (c) from the possible observation of 
a "magnetic field echo" in the auto-correlation function 
of conductance fluctuations. 

The latter prediction is made on the basis of the fol- 
lowing theoretical observation. After the unitary trans- 
formation (2), the leading spin-orbit term a*^ has pre- 
cisely the same spatial dependence as the vector potential 
due to the orbital magnetic fleld (in a symmetric gauge). 
Hence, can be regarded as an effective magnetic fleld 
perpendicular to the plane of the quantum well and of 
opposite sign for the two spin directions. A suitably cho- 
sen orbital magnetic fleld can balance this effective fleld, 
or change its sign for one spin direction, leading to a par- 
tial re-appearance of weak localization and to a "mag- 
netic fleld echo" in the auto-correlation function of con- 
ductance fluctuations. The magnetic fleld echo appears 
in the magnetoconductance autocorrelation function for 
a magnetic fleld difference 
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(8) 



The shift of weak localization peaks is by half this fleld 
strength. 

According to the Onsager's relations, the two-terminal 
conductance of a non- magnetic system is symmetric with 
respect to reversal of the magnetic field B. In the absence 
of spin-orbit coupling, the symmetry with respect to the 
reversal of the perpendicular magnetic fleld component 
B3 would persist even in the presence of Zeeman split- 
ting caused in a 2D electron gas by the in-plane magnetic 
field B|| (excluding effects arising from the finite width of 
the quantum well.^^) However, in the presence of even a 
weak SO coupling, the exact symmetry exists only if both 
B3 and B|| are reversed simultaneously, whereas the re- 
versal of the perpendicular component only generates an 
antisymmetric contribution to the magnetoconductance 
with variance 



var[G(B3,S||)-G(-S3,B||)] 



(iVA/27r)^ 



(9) 

where N is the number of channels in the leads, A is the 
mean level spacing in the dot, and {G) is the ensemble- 
averaged conductance of the quantum dot. Since the 
fluctuating antisymmetric contribution the conductance, 
G{B3 , B|| ) — G(— , S|| ) is linear in the in-plane magnetic 
fleld, £^ w gfigB^^, it dominates over the conductance 
asymmetry caused by inter-subband mixing due to the 
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FIG. 1: Illustration of the quantum dot with stub. Energy 
dependent scattering from the quantum dot in the presence of 
spin-orbit coupling and a magnetic field is modeled as scatter- 
ing from a chaotic quantum dot with a parameter independent 
scattering matrix U, with a stub described by a reflection ma- 
trix R which depends on energy, spin-orbit coupling and the 
magnetic field. 



in-plane magnetic field for intermediate magnetic fields 
strengths/'' which is proportional to Bf. 



III. AVERAGE CONDUCTANCE 

The starting point of our calculation is the Landauer 

formula for the two-terminal conductance G of the quan- 
tum dot at zero temperature, 



G 



2e2 N1N2 e 



—trSAS^A, 
h 



h N 

where the diagonal matrix A has elements 



" - 1 -N^/N, j 



Ni + 1,...,N. 



(10) 



(11) 



As a matrix of complex numbers, the scattering matrix 
S has dimension 2N = 2{Ni + N2); it may also be seen 
as an iV X TV matrix of quaternions, which are 2x2 
matrices with special rules for complex conjugation and 
transposition.^* 

In order to calculate the average of the conductance at 
Fermi energy e and magnetic field B, or the covariance 
at energies s and e' and magnetic fields B and B', it is 
sufficient to compute the average 



{Skl;nv{s, B)5fe';';^',y/(e', B')*), 



(12) 



where roman indices refer to the propagating channels 
in the leads, and Greek indices refer to spin. Hereto, 

we need a statistical description of the scattering ma- 
trix S in the presence of spin-orbit scattering and for 
arbitrary values of the magnetic field B and Fermi en- 
ergy e. Within the "random scattering matrix approach" 
of random-matrix theory, this is provided by the "stub 



model". -'^^ In this model, a fictitious "stub" is attached to 
the quantum dot, see Fig. 1. If the number of channels 
in the stub is M — N, the scattering matrix S can be 
written 



5 = pc/(i - g^i^Q^/)-lp^ 



(13) 



where U is the M x M scattering matrix of the quan- 
tum dot (without stub), R is the (M — iV)-diniensional 
reflection matrix of the stub, and P and Q are TV x M 
and (M — N) x M projection matrices, respectively, with 
Pij = Sijl, and Qij = 6i+N,j'^- (Here 1 is the 2x2 
unit matrix in spin space.) The area and width of the 
stub are chosen such that (i) the dwell time in the stub is 
much larger than the dwell time in the quantum dot and 
(ii) the time for ergodic exploration of the dot plus stub 
system is much shorter than the time for escape into one 
of the two leads. The first condition ensures that all de- 
pendence on Fermi energy, magnetic field, or spin-orbit 
scattering rate will be through the reflection matrix R, 
so that the matrix U can be taken at a fixed reference 
energy e = and at zero magnetic field and spin-orbit 
scattering rate. Since scattering from the quantum dot is 
chaotic, this implies that elements of U are proportional 
to the 2x2 unit matrix in spin space 1 and that U can 
be chosen from Dyson's circular orthogonal ensemble of 
random matrix theory. The second condition, which re- 
quires M ^ N, ensures that electrons explore the phase 
space of the combined dot plus stub system ergodically 
before they exit to the leads, so that the distribution of 
the scattering matrix S remains universal and is unaf- 
fected by the addition of the stub. Indeed, the paramet- 
ric dependence of the scattering matrix described by the 
stub model has been shown to be the same as that de- 
scribed by the Hamiltonian approach of random-matrix 
theory. 

We take the reflection matrix R of the form 



R = exp 



27ri 
MA 



{s-n') 



(14) 



where A is the mean level spacing of the closed dot and 
Ti' is an (M — A^) dimensional matrix which describes 
the effect of the magnetic fleld and spin-orbit scattering, 

A 

n = — [iX{xt + a^CTs) + ia {Aiai + A2(T2) 

-h-a + b^Bhas)]. (15) 

Here X, Ai and A2 are real antisymmetric matrices, with 
(trXX^) = M2 and (tr A^J) = while Bh is a 

real symmetric matrix with (^tr Bf^^ — M^. The dimen- 
sionless parameters x, a± , a, b, and b± correspond to 
the energy scales deflned in Eq. (6) as 



a\ = Tref/A, 

= Tref /A, 

h = TTS^/A, 

bl = nel/A. 



(16) 
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We define the corresponding quantities a;', 6', and b'^ 
when calculating a correlator between scattering matrix 
elements at different values B and B' of the magnetic 
field. For weak uniform spin-orbit scattering, one has 
the strong inequalities a <^ a±, b± ^ b. 

To calculate the correlator (12), we expand S in pow- 
ers of U using Eq. (13) and integrate U over the space 
of unitary symmetric matrices using the diagrammatic 
technique of Ref. 20. To leading order in 1/M and 1/A'' 
we can take the elements of U to be Gaussian random 
variables with mean zero and with variance (UijU^i) = 
{SikSji + 5ii6jk)- Then only ladder diagrams and 
maximally crossed diagrams contribute to the correla- 
tor Eq. (12). Summing the contributions from these two 
classes of diagrams yields 



where 



D ^ (Ml 1 - tr i? «) R'^) ^ 
C = {Mt®t-%xR®R'*)~^ . 




0.4 
0.2 



-0.2 
-0.4 
-0.6 



-4 

— I — 



(C) h=0, a=.5, aj^=3 



-I 1 1 — 

■ ■ ■ b=0 

— b=3 

— ■ b large 



(18) 



x/N 



1/2 



Here R' is given by Eq. (14), with x, b, b± and e re- 
placed by x', b', b'j_ and e', respectively. The super- 
script * denotes the quaternion complex conjugate of R', 
T = 1 (8) (72, and the trace is taken over the channel 
indices only. The tensor multiplication in Eq. (18) has 
a reverse-order multiplication for the Pauli matrices in 
second place. 



FIG. 2: Dependence of the quantum-interference correction to 

the conductance, SG, on the perpendicular magnetic field x. 
The conductance is measured in units of (2e^ / h)NiN2 / {Ni + 
N2Y' . Panel (a) shows 5G for ax = 3, a = 0, 6x = and 
two values of the parallel magnetic field, & = and 6 = 3, 
and for the limit 6 — > 00. Panel (b) shows 5G for a± = 3, 
6 = 6x = 0, and three values of the spin-orbit parameter 
a: a = 0, a = 0.5 and a — »■ 00. Panel (c) shows 5G in the 
crossover regime, with a± = 3.0, a = 0.5, and for the same 
values of the parallel magnetic field as in (a). In all three 



{<Ji (8> o-j){ai> (g) ay) = {aiai') (g) {cryaj). (19) Plots. we have set &x = 0. Choosing 6x < & is appropriate 



The two contributions C and D are the equivalents of 
cooperon and diffuson in the conventional diagrammatic 
perturbation theory.^^ Taking the limit M — > 00 and 
using Eq. (14) we find 



for small quantum dots with L <C Ai,2. A small nonzero value 
of 6x has only a slight effect on the shape of the curves, its 
general effect being to further reduce weak antilocalization. 



7V1 IH ^ -t (g) 1 



Nt^t^ ^-^ -t O 1 

+ ^a^(o-i(8)l-l(8ic7i)^ 

+ ^a^(o-2 01 -10(72)^ 

-I- ib • (<r (g) 1 - 1 O <t) 

1 



-(/ic73 (g) 1 - /I'l 00-3)^. 



(20) 



+ + a;')l ® 1 + "-lCo-s 1 - 1 (g 0-3)]^ 



.ai)2 



1 



a^(c72 (g 1 - 1 (gcr2)^ 



ib • (<r (g H- 1 (g <t) 



1 



(21) 



Setting a;' = a;, b = b' = bei and b'_^ = b±, we find the 
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average conductance, 

2e2 N1N2 



{G) 
SG 



h Ni+N2 



NiN2 



h N1 + N2 



+ dG, 



(22) 



N1N2 



h N1 + N2 V4a2 + Fc 



+ 



AGcb^ + G^Fc - l6alFcx^ 

where we abbreviated 

Fc = N + 2{x'' + bl) 

Gc = N + 2{x'' + a^+al). 



(23) 



(24) 



The expression for the average conductance simpUfies 
considerably when either the perpendicular magnetic 
field or the parallel magnetic field is zero. In the first 
case, b = b± = 0, one finds 



SG 



N1N2 
'hNi + N2 



1 



1 



N + 2x'^ 
1 



4a2 + AT + 2a;2 



^ 2a2 + iV • 



2{a±±xf 



(25) 



This simple result can be understood recalling that the 
spin-orbit coupling term a± acts as an efi'ective perpen- 
dicular field with opposite sign for up and down spins. 
When the applied perpendicular field x exactly cancels 
this effective field, i.e., when a; = ±a±, time-reversal sym- 
metry is "restored" for one spin direction, leading to a 
weak localization-like correction to the conductance cen- 
tered at a; = ±a±. A non-zero value of the spin-orbit 
coupling term a leads to a weak antilocalization peak at 
a; = 0, while a finite parallel magnetic fields causes a 
weak-localization peak at zero perpendicular field, as is 
illustrated in Fig. 2b and a, respectively. Both a parallel 
magnetic field and the spin-orbit coupling term a sup- 
press the features at x = ±a±. (The case a ^ a± is 
applicable when the spin-orbit coupling is not uniform, 
as discussed in the introduction.) Figure 2c shows an ex- 
ample for the dependence of the conductance G on the 
perpendicular field in the crossover regime when all pa- 
rameters (a_L, a, and b) are important. 

In the special case when the perpendicular magnetic 
field, x, is zero, we find 



SG 



N1N2 



1 



h N1+N2 [N + 2a^ + 2al 
2a^ 



2al 



+ 



462 + {N + 2h\){N + 2a? + 2a\) 
1 



N + Aa? + 2b\ 



In the case of uniform spin-orbit scattering, this may be 
further simplified using a?' <C and b\ <C 6,^ 



5G 



N1N2 



hNi+N2 
1 



2al 



+ 



■N 



(27) 



4a^ + 2b']_ + N 4b^ + + 2Nal 



At finite temperatures, dephasing will lead to a further 
suppression of weak localization. Dephasing can be in- 
cluded in the theory presented here by the introduction 
of a fictitious voltage probe. 22'23,24 rpj^j^ amounts to the 
replacement N ^ N + h/Tci,A in the final results (23)- 
(27), where is the dephasing time and A is the mean 
level spacing in the closed quantum dot.^^ 



IV. CONDUCTANCE FLUCTUATIONS 



Conductance fluctuations are described by the covari- 
ancc of the conductance at two different values e and e' 
of the Fermi energy and at two different magnetic fields 
B and B', 



cov[G(e,B),G(£',B')] 



{G{e,B),G{e',B')) 
-{G{e,B)){G{e',B')). 



(28) 



For this calculation, we need to know the average of a 
product of four scattering matrix elements. However, 
if the conductance is expressed in terms of the scatter- 
ing matrix S using Eq. (10), to leading order in 1/N, 
the scattering matrix elements may be considered Gaus- 
sian random numbers, and the average of four scattering 
matrix elements can be factorizcd into products of pair 
averages of the form of Eq. (12). 2^ One thus obtains 



cov [G(e,B),G(e',B')] 



N1N2 



Vc). 



No 



(29) 



(26) where 
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Vd = ^ D ni D 1,1^1.1,^ (30) 

1(62 - 6'2)2 + 2ifD(-4a26&' + (&2 + 6'2)L^) + (l2^ _ ^a^){Kl - Aa\{x - x'f)f ' 
Vc = Yl C^C^), (^C-^),'.';,. (31) 

_ 2c 

1(62 - 6'2)2 + 2Kci^a^bb' + (62 + 6'2)Lc) + (i^ - 4a4)(ii:2 _ 4al{x + x')2)|^ ' 

Here we abbreviated 

Ed = 2 iKDi'la" - 1%) + 4a%b' - ^^(6^ + 6'2)|' + 2 17^0(6^ + 6'^) + kILd - Aa^Loix - x'f\^ 

+ 8ai(a; - x'f ((6 + b'f + |Ld + 2a2f ) ((6 - 6')^ + |Ld - 20^1^) + 8 |i<:i566' + a2(i<:|, - Aa\{x - x'f)\^ 
+ 4 |6(62 - 6'2) - KD{2a%' - 6Ld)|^ + 4 |6'(6'2 - 6^) - Koi^aH - b'LD)\^ + 8 |a'(62 + 6'^) - 66'Ld|' 

Sc = 2 |i^c(4a4 - L^) - 40^66' - Lc(&' + «>")r + 2 |irc(&' + &") + K^Lc - AalLc{x + x'f\^ 

+ 8oi(a; + x')' ((6 - b'f + \Lc + 2a2f ) ((6 + b'f + \Lc - 2a^\^^ + 8 \Kcbb' + a^i-K^ + 4al{x + x'f)\'^ 
+ 4 |6(62 - 6'2) + ii:c(2a26' + 6Xc)f + 4 |6'(6'2 - 62) + ii:c(2a26 + 6'Lc)f + 8 |a2(62 + 6'2) + 66'Lcf 

and 





= N + 


2m{e - 
A 


e') 


+ 


l{x-x'f + 2{a'' + al) + \{b^ + b'^f 


Ld 


= N + 


2m{£ - 
A 


e') 


+ 


i(x-x')' + 2a2 + 1(6^-61)% 


Kc 


= N + 


2m{£ - 
A 


e') 


+ 


^-{x + x'f + 2{a'' + al) + \{b^-b'^f 


Lc 


= N + 


2m{£ - 
A 




+ 


i(x + x')' + 2a2 + 1(6^ + 61)'. 



A. Variance 



The variance of the conductance at zero temperature is obtained from Eq. (29) by setting e' = e and B' = B, 



var G 



NlNl 



h? {Ni + N2f 



' ^ J_ 862 + G2, + (4a2 _^ 1 
.7V2 + ^ {Aa?GD + 462 + GoNf ^ {Aa? + Fcf 

2Fc 



+ 



Gc(462 + GcFc) - lQa\Fcx-^ 
where Fc and Gc are given by Eq. (24), and 

Gz3 = iV + 2 (a^ + ai + 6^) . 



(462 + GI + IQalx^y + 640^x2(^2 , d 
{Gci^b^ + GcFc) - l&aiFcx-^ f 



(32) 



(33) 



Simphfications occur in the hmits of zero parallel or perpendicular field, and for large perpendicular field. In the 
absence of a parallel field, 6 = 6_l = 0, one finds the variance 



vaxG 



/l2 (iVi + 7V2)2 
+ 



1 



+ 



1 



+ 



1 



(2a2 + AT + 2(a_L - a;)2)2 (iV + 2a;2)2 (4a2 + TV + 2a;2)2 

11 2 " 



(2a2+iV + 2(ax+a;)2)2 iV2 (4a2 + iV)2 (2(a2 + a^) + iV)2 



(34) 
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while in the absence of a perpendicular magnetic field, x = 0, one has 



varG = 



{Ni+N2y [N^ (4a2 + 262_ + iV)2 [2{a^ + al + bl) + N]"^ 4^2 + (4a2 + 7V)[2(o2 + a^) + iV] 

4(a2 + a^)2 4(a2 _ a^)2 



+ 



462 + 7V(2(a2 + a^) + A^) (4^2 + N[2{a^ + a^) + N]f (462 + (4a2 + iV)[2(a2 + a^) + N]f 



.(35) 



In the last equality we also used b± <C 6^. In the limit of large perpendicular field, x ^ 1, one has 



var G 



/i2 (TVi + N^)^ 



1 



1 2 
iV2 ^ 462 + (4o2 + N) [2(a2 + + ^2 ) + ^v] 

A{-a^ + al+blf 



[2{c 



il + 61) + NY (462 + (4a2 + N) [2{a^ + al + 6i) + N])' 



(36) 



When both the perpendicular and the parallel field are large, i.e. x^ ,b ^ a\, ,b\, N a particularly simple expression 
for the variance results 



var G 



{Ni + N2f 



1 



1 



(iV + 2a2 + 2ai + 26i)'_ 



5 
4 

O 3: 



^ 2 







! — — 1 — 


— 1 — ' — 1 — ' — 1 — ' — ; 

— b=0 

- ■ b=3 - 
■ ■ ■ b large '. 


\ ^ 




:a 

* 




■ . 1 . 






1 



2 

x/N 



3 

1/2 



FIG. 3: Conductance fluctuations varG, measured in units of 

as a function of the perpendicular 
magnetic field x, for a_L = 3, a = 0, &x = 0. From bottom to 
top, the curves are for 6 = 0, 6 = 3, and for the limit 6 — > 00. 



direction is real, not complex. A parallel field and the 
spin-orbit coupling term a suppress this feature (data for 
a not shown). Both effects also reduce the variance at 
zero perpendicular field by a factor of two. 

At finite temperatures, both the effects of thermal 
smearing and dephasing need to be taken into account 
for a calculation of the variance of the conductance. De- 
phasing (with dephasing time r<^) is incorporated in the 
results (29)-(37) by the substitution N ^ N + h/r^A. 
The effect of thermal smearing requires an integration 
over the energies e and e', 



varG(T) 



f, , ,df df 



coY[G{e),G{e')], (37) 



where / is the Fermi distribution function. The integra- 
tions over energy needs to be done by numerical methods 
in all but a few special cases. We refer to Rcf. 26 for 
details. Equation (37) was used for a quantitative com- 
parison of theory and experiment by Zumbiihl et al}^ 



Symmetry with respect to perpendicular 
magnetic field inversion 



The variance of the conductance in the presence of 
spin-orbit coupling is shown in Fig. 3, for three different 
parallel field strengths. An increase of the conductance 

fluctuations is observed around x = ±a^. This feature 
is explained noting that at a; = ±a_L the perpendicular 
magnetic field exactly cancels the effective magnetic field 
due to the spin-orbit coupling term a_\_ for one spin di- 
rection, so that the effective Hamiltonian for that spin 



The presence of spin-orbit coupling has interesting con- 
sequences for the symmetry of the conductance under re- 
versal of the perpendicular field. Although, in general, 
the conductance does not change when the total mag- 
netic field is reversed, 6 —6, 6^ ^b±, and x — x, 
the conductance need not be symmetric under reversal of 
the perpendicular component x only. In order to study 
the symmetry of the conductance under reversal of the 
perpendicular component x of the magnetic field, we con- 
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sider the correlator 



Q{x) = 



(G(x)G(-x)) - (G(x))^ 
var G{x) 



(38) 



for ^ N. Perfect symmetry under reversal of the per- 
pendicular component of the magnetic field only implies 
Q{x) = 1, while the absence of any correlations between 
G{x) and G{—x) implies Q{x) = 0. 

Note that the invariance of the two-terminal conduc- 
tance conductance under a complete reversal of the mag- 
netic field implies that 

{G{x, b)G{-x, b)) = {G{x, b)G{x, -b)) , 

so that Q can also be defined as a correlator for conduc- 
tances at opposite values of the parallel magnetic field, 
keeping the perpendicular field constant. 

In the absence of an in-plane magnetic field, reversal of 
the perpendicular magnetic field component is the same 
as reversal of the total magnetic field, so that G is triv- 
ially symmetric in x. In the absence of spin-orbit scat- 
tering, the sole effect of the parallel magnetic field is to 
lift spin degeneracy, and symmetry of G follows from the 
fact that the conductance for each spin species is sym- 
metric in X. In the presence of both spin-orbit scattering 
and an in-plane field, there is no general symmetry that 
requires that G is symmetric under reversal —B^, 
however the measure of observable asymmetry depends 
upon to what cross-over symmetry class the generated 
random matrix theory ensemble belongs. 

Indeed, calculation of the correlator Q{x) using our 
general result (29) for the covariance of the conductance 
shows that the symmetry of the conductance with re- 
spect to reversal of the perpendicular component of the 
magnetic field is destroyed when either 6 7^ and a ^ 0, 
b± and a± 0, ot b± and a 7^ 0. In the case 
a = and b± = 0, the conductance remains symmet- 
ric in X, even if a_L 7^ and b ^ 0. In that case, the 
effective Hamiltonian (4) still has the extra symmetry 
(JxHax = H, which ensures the symmetry with respect 
to the perpendicular magnetic field reversal.^ Symmetry 
in the perpendicular component of the magnetic field is 
also preserved if the parallel magnetic field is so large 
that both b, b±^ ^ 0, but < and <C N, see Ref. 
7 for a symmetry argument. Figure 4 shows plots of Q{x) 
for the two minimal cases o = or 6_l = 0. 

The general parameter dependence of Q{x) for x'^ ^ N 
can be obtained from Eq. (29). Simple expressions are 
found in the cases b ^ N and a\ ^ N, for which we 
have 



Q 



(4a2 



2bl 



Ny 



and for 6 > iV, while still 6_l = 0, 
7V2 



Q = 1- 



(4a2 + Ny 



(39) 



(40) 



2[(a2 + ai)2 + (a2 + ai + Ar)2] 



1 

0.8 
00.6 
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0.2 

0, 



1 ' 
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FIG. 4: The correlator Q{x) = {{G{x)G{-x)) - 

(G(a;))^)/varG'(a;) for x'^ > iV. In panel (a), Q is shown 
as a function of parallel field b, with = = 0. In panel 
(b), Q is shown as a function of spin-orbit coupling ax, with 
6 = 5 and a = 0. 



On the other hand, iib ^ N and b± ^ N, the correlation 
function Q is given by 



Q = 



(2a2 + 2al + Ny 



(41) 



confirming that G remains a symmetric function of x if 
b and b± arc large but a and a± are not. 

Experimentally, the component of the conductance 
fluctuations that is antisymmetric with respect to the 
reversal B3 —B^ in the presence of in-plane magnetic 
field can be used as a test to the strength of spin-orbit 
coupling, even if 6 <C iV. This is because an alterna- 
tive mechanism of time-reversal symmetry breaking by 
an in-plane magnetic field via the inter-subband mixing^^ 
generates anti-symmetric condiictance fliictuations in the 
third order in the field only, thus 1 — Q{x ^ N) ^ B^^. 
In contrast, the spin-orbit-coupling induced asymmetry, 
which can be analyzed by expanding Q{x) with respect 
to 52 and 6^, has 
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Q = l-4[4a^b^i32a*ai + 24a^a'iN + 6a^^N^+4alN^ + N^) 

X [Ar(4a2 + iV)(2ai + N){lQa^a\ + ia'a^N + 2a\N^ + 2a\N'^ + N^)\~^ , (42) 



for the case of uniform spin-orbit coupling. In the limit 
of a dot with a large number of channels N a\ but 
still in the universal regime A'' ^ i^Th/A, this simplifies 
to 



Q = 1 - l6a%^N- 



(43) 



Replacing the dimensionless spin-orbit coupling and mag- 
netic field a±_ and 6^ by their dimcnsionful counterparts, 
one arrives at the result (9) advertised in the introduc- 
tion. 



C. Autocorrelation function and the "correlation 
echo" 

Correlations between the conductance G at different 
values of the perpendicular magnetic field x are studied 
through the conductance autocorrelation function 



x{x,x') = {G{x)G{x'))-{G{x)){G{x')). 



(44) 



For large perpendicular fields, x'^,x''^ 3> N, the autocor- 
relator x{x,x') depends on the difference {x — x') only. 
Since time-reversal symmetry is fully broken for such 
large magnetic fields, any features in the autocorrelator 
are signatures of the spin-orbit coupling. 

For zero parallel field, 6 = 6_l = 0, one finds from Eq. 
(29), 



X(0 = 



2e2 N1N2 



+ 



h N1 + N2 
1 



+ 



{2N + 8a2 + $2)2 
1 



(27V -h $2)2 (27V 4a2 (-2a_L 0^)^ 
1 



(2iV + 4a2 + (2ax+0')' 



(45) 



where = x — x' . 

A plot of x(x — x') for different values of a± is shown 
in Fig. 5. Remarkably, the conductance autocorrelator 
shows an "echo" at |a; — = 2a ±. This echo appears 
even for large values of a±, for which the conductance 
correlations have decayed for intermediate values of the 
magnetic-field difference, see Fig. 5a. The magnetic-field 
echo has the same origin as the reappearance of the weak 
localization and conductance fluctuations at flnite val- 
ues of the perpendicular magnetic field x: The spin-orbit 
term a± acts as a homogeneous perpendicular field with 
a different direction for the up and down spins (defined 
with respect to the axis 63 perpendicular to the plane of 



X 




(x-x')/N 



FIG. 5: The correlator x{x ^ &s a function of the dimen- 
sionless perpendicular magnetic field difference {x — x'), for 
three values of the spin-orbit term a±: a± = 0, a± = 1, and 
a± = 3. Curves in panel (a) are for a = and b = b± = 0. 
Curves in panel (b) are for o = 0.5 and b = b± = 0. 
Panel (c) shows the effect of the parallel field, 6 = 2, with 
a = b± — 0. In all three panels, x{x — x') is measured in units 
of [{e'/h)NlNi/{Ni + N2f]^. 



the quantum well). In the presence of an external per- 
pendicular magnetic field x, Up spins and down spins 
experience effective perpendicular magnetic fields x + a± 
and X — a±_ , respectively. If the magnetic field difference 
X — x' equals ±2a^, one spin species at field x moves in 
the same effective perpendicular field as the other spin 
species at field x' . Hence, in the absence of a parallel field 
and a spin-orbit term a, the conductance has identical 
fluctuating contributions at x and x' . However, only one 
spin species is involved in the "echo" — the autocorrela- 
tion function reaches only one half of the total variance 
of the conductance; the other spin species have uncor- 
related contributions to the conductance at the external 
flelds x and x' . A magnetic field in the plane of the quan- 
tum well and the spin-orbit coupling term a suppress the 
echo, since they cause scattering between the ±63 spin 
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directions, see Figs. 5b and c. 



V. EFFECT OF A NON-UNIFORM SO 
COUPLING. 

The expressions of the previous two sections are valid 
for arbitrary magnitudes of the dimensionless spin-orbit 
rates a± and a. In the case of uniform spin-orbit cou- 
pUng, one has the strong inequahty a <C a_L. In general, 
the spin-orbit coupling may be non-uniform throughout 
the quantum dot, e.g., result of a spatial modula- 
tion of the asymmetry of the quantum well or as a result 
of a modulation of the electron density (and therefore 
confining potential width) across the dot. The main ef- 
fect of a non-uniformity of the spin-orbit coupling is an 
increase of a,^^ thus promoting the crossover into the 
standard symplectic ensemble. In this section, we'll give 
a quantitative estimate of this effect for a special choice 
of the non-uniformity. Another effect of a non-uniform 
spin-orbit coupling is a rescaling of a±, and will not be 
considered here. 

Assuming that the non-uniformity of spin-orbit cou- 
pling affects only its strength without changing the 
generic form of the Hamiltonian, we describe its effect via 
coordinate dependent spin-orbit coupling lengths Aj~2(r), 



1 
2m 

0-2 



n = — (p - eA 



2Ai(r) 2A2(r) 



62- 



(46) 



In the case that 



/3 = V X a = y9,,Ar^(r) + ^d,^X^\r) ^ 0, 

it is not possible to eliminate the spin-orbit coupling from 
the Hamiltonian in the first order in spin-orbit coupling 
constant by a unitary transformation. Qualitatively, this 
explains why the spatial variation of the spin-orbit cou- 
pling length gives rise of a term with the symmetry of 
a|| . For a quantitative estimate of how much the spatial 
variations of Ai,2 contribute to the energy scale £||, we 
note that it is still possible to find a transformation U 
such that the resulting non-Abelian vector potential in a 
locally rotated spin frame. 



CTl 



a = - i[/tVJ7 - U^faU, 
would satisiy the conditions 

V-a = 



(47) 



(48) 



in the interior of the quantum dot and n • a = 0, on 

the dot boundary, where n is the unit vector normal to 
the dot's boundary, at least, in the lowest order in SO 
coupling. Such conditions would enable one to satisfy 
boundary conditions for Cooperon and diffuson propa- 
gators from diagrammatic perturbation theory, which is 



needed in order to make the zero-dimensional approxi- 
mation necessary for the use of random matrix theory. 
This program is carried out in the appendix for a circu- 
lar quantum dot with a special choice of the coordinate 
dependence of the spin-orbit coupling lengths Ai^2- 

Here, as an example, we consider here the case where 
the spin-orbit coupling is stronger in the center of a quan- 
tum dot (where density is higher and the well is narrower) 
than at the edges (where density is lower and the well is 
broader), 



1 

Al,2 



1 

A^ 



Af 



(49) 



where L is the typical size of the dot, without making a 
restriction of the size of the dot. We perform a unitary 
transformation 



U = exp 



.Xia2 — X2(J\ 

' 2A^ 



3At VZJ 



(50) 



which rotates to a local spin reference axis. With this 
transformation, the matrix vector a becomes, to lowest 
order in the inverse spin-orbit coupling length, 



a = -z[/tvC/-a-FO(A-^) 
= SA^I"^^'-]- 



(51) 



We conclude that (a) the spin relaxation and the 
crossover into symplectic ensemble start, indeed, in the 
lowest order in the non-uniform part of the spin-orbit 
coupling and (b) for the non- uniform coupling (49), the 
functional form of the perturbation after transformation 
is equal to that of the term ap in Eq. (4). The latter 
observation allows us to express the contribution of the 
non- uniformity to the energy scale £11 as 



e|| = ;^K'Ki^Th 



4Af2 



(52) 



where the geometry-dependent coefficients k and n' are 
the same as those defined for the spatially homogeneous 
spin-orbit coupling in Eq. (6). 



VI. CONCLUSIONS 

In this paper we calculated the average and the fluc- 
tuations of the conductance of a chaotic quantum dot 
in the presence of Bychkov-Rashba and (linear) Dressel- 
haus spin-orbit coupling and a magnetic field and stud- 
ied the symmetry of the conductance under reversal of 
the perpendicular component of the magnetic field. The 
calculations were done to leading order in 1/A'', where 
N = Ni + N2 is the total number of open channels in 
the leads connecting the quantum dot to the reservoirs. 
After a unitary transformation, the spin-orbit term in 
the Hamiltonian and the Zeeman coupling to the exter- 
nal magnetic field are transformed into four terms that 
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all have different symmetries/ The corresponding energy 
scales are the conventional Zceman energy — HsgB 
and three energy scales that depend on the spin-orbit 
coupling: two spin-orbit scales ej^ and and the sec- 
ond Zeeman energy scale e"^. With the help of the theory 
presented here, the energy scales e^°, e^^°, and e'j_ can be 
obtained from a measurement of the conductance auto- 
correlation function of the quantum dot. 

Our results provide several routes to a direct measure- 
ment of the product of the spin-orbit coupling lengths 
Ai and A2. The product A1A2 can be foimd from a 
"magnetic-field echo" in the conductance fluctuations 
as a function of the perpendicular magnetic field. The 
magnetic-field echo appears at the magnetic-field differ- 
ence 



AB = 



eAiAi 



(53) 



The same magnetic-field scale splits the weak-localization 

peak in the absence of a parallel magnetic field. The 
ratio A1/A2 can be found from a measurement of the 
energy scale in a roughly circular quantum dot. In this 
case, (A2/Ai)^ is equal to the ratio of the energy scales 
e"^ for the in-plane magnetic field in the crystallographic 
directions ei = [110] and 62 = [110]. Measurement of e||° 
is less suitable to determine A1/A2, as it may be affected 
by non-imiformities in the spin-orbit scattering rate. 

The theory presented here has been used for the inter- 
pretation of measurements of the conductance of quan- 
tum dots in GaAs/GaAlAs hetcrostructurcs with spin- 
orbit coupling as a function of parallel and perpendicu- 
lar components of the magnetic field by Zumbiihl et al}^^ 
The quantum dots used in the experiments of Ref. 15 
were built of the same material, but they varied in size, 
ranging from L/X = 0.27 to 0.64. Here A = (AiA2)^/^ is 
the geometric average of the spin-orbit lengths Ai and A2. 
Whereas the conductance of smaller dots had a minimum 
at zero perpendicular magnetic field, the conductance of 
larger dots had a (local) maximum at zero magnetic field. 
Prom a fit to our expression for the average conductance 
at zero parallel field, Eq. (23), values for the experimen- 
tally unknown parameters A = 4.4/xm, the decoherence 
time T0 and a geometrical factor, see Eqs. (6), were ex- 
tracted. Measurement of the average conductance in the 
presence of a magnetic field parallel to the plane of the 
two-dimensional electron gas was in good agreement with 
our Eq. (23), without additional fit parameters and up 
to a parallel field of about 0.3 T.^^ Upon inclusion of the 
effects of time-reversal symmetry breaking by a strong 
parallel field due to the finite extent of the electron wave- 
function across the heterostructure,^'' the authors of Ref. 
15 could extend the agreement between theory and ex- 
periment is extended to even higher magnetic fields. 
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APPENDIX A: GEOMETRY-DEPENDENT 
COEFFICIENTS FOR A CIRCULAR QUANTUM 
DOT 

In this appendix, we calculate 

xG%,^{v„v,){E',B')) (AI) 

for a disordered circular quantum dot with spin-orbit 
scattering. Here and are the retarded and ad- 
vanced Green functions (inverses of e — W for £ just above 
and just below the real axis, respectively), respectively, 
and the indices /xi, /U2, vi, and 1/2 refer to the electron 
spin. Comparison of our result with the random ma- 
trix theory of Sec. Ill allows us to find the geometry- 
dependent coefficients and k" of Eqs. (6). 

In a closed dot of volume V and with diffusive dynam- 
ics, the diffuson ■D(ri,r2) obeys the equation 

277 

{ihw + n)l?(ri , V2) = y^S{ri - r2)l 1, (A2) 

where fko = e' —e and the differential operator H is given 
by 

H = {D/h){~ihVr,+a(E)i-l'E)a'f 

i/iz «) 1 - O /I'z- (A3) 

Here D = vpl/'^. is the diffusion coefficient in the quan- 
tum dot, the matrix vector potentials a and a' are the 
sum of a spin-independent contribution from the per- 
pendicular magnetic field and a spin-dependent contri- 
bution from the Bychkov-Rashba and Dresselhaus spin- 
orbit coupling, 



eBsy 


h 


2c 






h 


2c 


2X1 



^2 ei - 



0-2 ei 



/ eB^x h 
\^ ~ 2X2 



eB'^x h 
~2c 2A^ 



CTl 62, 



cri 62, 



(A4) 



while hz and /i^ represent the Zeeman coupling to the 
magnetic field, 
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Equation (A2) is supplemented with the boundary con- 
dition 

n ■ {-ifiVri + a ® 1 - a ® a')I?(ri , r2) = (A5) 

at the sample boundary, where n is the unit vector nor- 
mal to the boundary. 

To find X>, we shift to the basis of eigenfunctions of 
the diffusion equation (i.e., eigenfunctions of 11 in the ab- 
sence of spin-orbit scattering and a magnetic field) and 
calculate the matrix elements of 11 in that basis. For this 
procedure to work, it is important that the matrix vec- 
tor potentials a and a' do not change the Von Neumann 
boundary conditions (A5), i.e., that 



n • a = n • a 







(A6) 



at the boundary of the quantum dot. This requirement 
ensures that the eigenfunctions of the diffusion equation 
with Von Neumann boundary conditions can be used to 
construct the diffuson V. In order to fix the gauge, we 
also require that 



V • a = V • a' 







(A7) 



in the interior of the dot (London gauge). 

As discussed in Rcf. 7 and in the introduction of this 
paper, the requirements (A6) and (A7) are realized by 
performing a suitable unitary transformation U to the 
Hamiltonian, H U^HU. For Eq. (A2), such a unitary 
transformation implies 



n {u'^ ®u)T\{u ®u'^). 



(A8) 
(A9) 



We'll now perform an explicit calculation for the case 
of a circular quantum dot of radius L, the origin of our 
coordinate system being chosen at the center of the dot. 
In polar coordinates xi = r cos (p, X2 = r sin (p, the tensor 
eigenfunctions of the diffusion equation are of the form 

fnoir^^)) = CnoJo{rXno/L)(Ji® <7j, 

{rXnm/L) cos(m^) a-j (g) a^, 
{rXnm/L,) sm{m<j)) <Ji iTj(AlO) 



where i,j = 0,1,2,3, m = 1,2,..., Xr, 



1,2,..., 



is the nth positive root of Jm{x) = 0, and we used the 



notation ao = 1. The corresponding eigenvalues are 

Xnm = Dh{Xnm/Lf (All) 

and the normalization constants are 

Cnm = ^'^L'^{Jm{Xnrnf - Jrn+l{Xnrnf){^ + 5rno)- (A12) 

The smallest eigenvalue is found for n = m = 0, Aqo = 0; 
all other eigenvalues are larger than the Thouless energy 
Et^ = Dh/L^. 

For a circular quantum dot, the magnetic-field contri- 
bution to the matrix vector potential of Eq. (A3) already 
satisfies the requirements (A6) and (A7). For the spin- 
orbit contribution we use the transformation 



U 



exp 



.Xi(J2 .X2CTI 
I : l- 



2Ai 



. {x 



48A1A2 



2A2 
3L^) f x2ai 
\ Ai 



(A13) 



A2 



which fulfills the requirements (A6) and (A7) to leading 
order in L/ A. [Note: although the transformation of Eq. 
(2) is sufficient to satisfy the requirements to order L/\, 
the remaining term a|| of Eq. (4) does not obey Eq. (A7).] 
With this transformation, the matrix vector potential 
a and the Zeeman term hz change to 



ei 



+ 



683X2 hX20-3 



2c 

h 



4A1A2 

/ 2a;ia:2cri 
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+ 



16A1A2 V 

eBsXi 



+ 



Ai 
hxias 



+ 



{3x1 



.T?-L>2 



2c ■ 4A1A2 

h f ^X\X2(J2 



I6A1A2 V A2 



+ 



A2 



{7,x\ +x\- L'^)ai 



1 „ 1 



BiXi 
Ai 



+ 



Ai 

B2X2 
A2 



(A14) 



with similar changes for a' and /ig- 

As long as all relevant energy scales are much smaller 
than i^Th, it is sufficient to know the action of 11 on the 16 
eigenfunctions at /qq at n = m = 0. In tensor notation, 
this action is 



n 



Dh 
21? 

DHL^ 



/ e 

— ($ - $')! ® 1 + 7To(^3 O 1 - 1 O C73 

\ cfi 4A^ 



+ 



+ 1 



96A4 



CTl 



)CTl 



(^' 



2Ai 
1-1 



DHL^ / (T2 1 - 1 C72 



96A4 

,2 „2 r4 



2A2 



+ -| ) (1 CT3 - C73 1)' 
Al A2 



(A15) 
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where $ = B^nL"^ is the magnetic flux through the quan- 
tum dot and 



With this transformation, the matrix vector a is mapped 
to 



1 



2 ^Ini^ln 



1) 



0.0182 



is a numerical constant. Confining ourselves to the case 
A = Ai = A2 and comparing this result to Eq. (20), we 
can make the identifications 



b = 



4A2 

4AV 3L2A' 



hvpln 

UK' 



b± 



1.83 



(A16) 

(A17) 

(A18) 
(A19) 
(A20) 



For the geometry-dependent constants and k" in 

Eq. (6), this implies k = 2, k' = 1/3 and k" w 0.292. 

If the spin-orbit coupling is not uniform throughout the 
quantum dot, see Sec. V, the above estimates for a_L and 
a and the corresponding energy scales and ej'" need to 
be revisited. Here we consider the example of Eq. (49), 
^12^^) = K-^ + and construct the unitary 

transformation U that ensures that the conditions (A6) 
and (A7) are satisfied to lowest order in 1/A/, 



ei 



eB^X2 hX2'J3 

^iA^Af 
+ 



2c 

h 

8l2 



2xiX2cri , {3x2 



■xl 



L^)<T2 
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eB^xi , /;,.Ti(T3 



2c 



'2X1X2(72 {ix\ 

8L2 V \\ ^ 



xl 



(A22) 



Comparing Eqs. (A22) and (A14), we immediately con- 
clude that with a non-uniform spin-orbit scattering 
strength of the form (49) with A* = A^ = A2 one should 
make the identification 



(4(Af)2 ; 12L2A" 



(A23) 



In terms of the energy scale £m°, this implies 



U 



exp 



.xia2 .X2cri 
i-TTT «- 



2\\ 



2\% 

„2 



-I- ix\a2 



2 

{x'i +xl+ L^) 



IX20\- 



{x\ + xl + L2) 



8L2Af 



(A21) 



"II - 6^^^ V4(A02 I ' 



in agreement with the general relation (52). 
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